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Abstract 

A viscous-inviscid interactive procedure for 
subsonic flow is developed and applied to an axial 
coropresscr stage. Calculations are carried out on 
a two-dimensional blade-to-b!ade region of constant 
radius assumed to occupy a mid-span location. Hub 
and tip effects are neglected. The Euler Equations 
are solved by MacCormack's method, a viscous march- 
ing procedure is used in the boundary layers and 
wake, and an iterative interaction scheme is con- 
structed that matches them in a way that Incorporates 
Information related to momentum and enthalpy thick- 
nesses as well the displacement thickness. The 
calculations are quasi- three-dimensional in the 
sense that the boundary layer and wake solutions 
allow for the presence of spanwlse (radial) veloci- 
ties . 


Nomenclature 


B(x) 

Lower boundary of the cascade 'tolutlon 
region (Fig. 1) 

c 

Speed of sound 

C 

Blade chord length 

c ,c 
P V 

Specific heats 

e 

Total energy (Eq. 4) 

f 

Composite solution vector defined by Eq . 

c 

27 

fb(0) 

Boundary layer solution vector (Fig. 4) 

f.g 

Vectors defined in Eq. 23 


Vectors defined in Eq . 22 

F,G,H 

Vectors defined by Eq . 4 

(«b(°»o 

Vector obtained by evaluating g using 
values obtained from the boundary layer 
solution, fj^(O), (Eq. 27) 

G1,G2,G3, 

G4 

If 1 U if 

Components of vector G 

m £f 


41 2 3 , Components of vector 


h 

Enthalpy 

i.J 

Grid indices (Fig. 2) 

k 

Thermal conductivity 

P 

Pressure 

r,9.Z 

Cylindrical coordinates used in the vis- 


cous solution (Fig. 3) 

R 

Gas constant (Eq. 4) 

t 

Time 

T 

Temperature 

U 

Vector defined by Eq. 4 

u ,u 

Velocity components used in the invlscid 

s n 

solution (Fig. 1) 
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u,v,w 


x.y 

X,Y,Z 


Of, 3 
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s 


e 


5.C.11 


p 

♦ 


0 


Velocity components in the x and y 
directions (Eqs. 2,3) 

Velocity components in the and 
directions 

Velocity component in the and 

7] directions (Fig 3) 

Cartesian coordii ates used in the 
invlscid solution (Fig. 1) 

Cartesian coordinates used in Sec- 
tion 4 (Fig. 4) 

Cartesian coordinates used in the 
viscous solution (Fig. 3) 

Cascade coordinates (Fig. 1) 

Boundary layer thickness (Section 

4) 

Boundary layer thickness, nondimen- 
sionallzed with respect to C(Sec- 
tlon 3) 

Displacement thickness ^ 

Small parameter that is ^’(«) 
angle between B(x) and the x-dlrcc- 
tlon in the cascade solution region 
(Fig. 1) 

Viscosities 

Curvilinear coordinates used in the 
viscous solution (Fig. 3, Eqs. 11, 
12,13) 

Density 

Angle between the ^'Coordinate line 
and a cylindrical generator (Fig. 3) 
Angular velocity of the blades (rad/ 
sec) 


Subscrlptf 

REF 

WALL 

o 

6 


Denotes reference quantity 
Denotes evaluation at wall 
Denotes evaluation at y“0 
Denotes evaluation at >*6 
Denotes evaluation at infinity 
(Section 3) 


S uperscripts 


Denotes derivative (Eqs. 8,9) 


I . Introduc tion 

The flow in the blade passages of an axial 
compressor is quite complicated. In general, the 
flow is compressible, viscous, unsteady with re- 
spect to the blades, turbulent, and highly three 
dimensional. Furthermore, the flow may be either 
entirely subsonic or at least partially supersonic. 
Boundary layer separation may occur at several 
locations in a compressor blade passage. In addi- 
tion, any computational attempt to deal reai.'i'i- 
cally with such a flow will encounter these dxifi- 
culties in a region which i« geometrically complex. 

For these reasons, it is unlikely that a com- 
pletely realistic solution of the flow through an 
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entire axial compressor blade passage conid be un- 
dertaken In the near future. However, a great deal 
has already been accomplished by examining problems 
of a simplified nature. The numerical compressor 
calculations which have been undertaken have tended 
to examine certain aspects of the flow in a blade 
passage while ignoring other effects, even though 
the effects which are omitted frixn the calculation 
may be quite Important In practice. In this manner, 
it has been possible to reduce the original problem 
to one which is ir.athematically tractable, and to 
gain ome insight into the nature of such a flow. 

There have been several approaches to simplifi- 
cation. One approach, presented in detail by Wu^ , 
consists of the specification of a stream surface 
within the blade passage and the subsequent solution 
of the Euler equations on that stream surface. The 
surface can be specified either as an annular surface 
(blade to blade) or a meridional surface (hub to 
shroud), and by confining the calculation to a two- 
dimensional region in this fashion, it is possible 
to introduce the passage geometry into the calcu- 
lation while retaining the numerical benefits of a 
scalar stream function. This technique has been 
used by Katsanls? and Katsanis and McNally^, among 
others. The computer codes of references 2 and 3 
are well established, and are currently used in 
compressor design. 

Boundary layer calculations have been done on 
compressor blades and on passage cndwalls^*^. Simi- 
lar to boundary layer calculations are the viscous 
nutrching procedures, which are used to solve para- 
bolized Navier-S tokes equations^*®. Viscous march- 
ing procedures are currently being applied to flow 
in turbomachinery. 

Another popular simplification of compressor 
flow is its idealization as flow through a cascade 
of airfoils. Many different nimerical solutions 
have been carried out in cascade geometries. Per- 
haps the most studied set of equations with reference 
to cascade geixnetrles are the Euler equations, and 
a popular approach to their numerical solution has 
been through time marching techniques . Hiese 
techniques, of which MacCormack's method^^ has been 
the most widely used, owe their popularity to several 
factors. They are computationally efficient, they 
can be used for both subsonic and supersonic flow, 
and they are not subject to some of the limitations 
of simpler solution methods, such as irrotat ional ity 
and two-dimensional ity. And certain recently devc- 
loped time marching a Igor i thms , which are 
applicable ... the solution of the Ei ler equations 
and Navler-S tokes equations, appear to be quite 
promising for increased computational efficiency. 

It is likely that these new algorithms, or variants 
of them, will be used in the near future to carry 
out compressor calculations of increasing sophisti- 
cation. 

While the preceding survey of numerical compres- 
sor calculations is by no mi>ans complete, it servos 
to demonstrate the diversity of approaches to the 
overall problem, which is too difficult to be 
attacked in a more straightforward manner. The 
present investigation is primarily concerned with 
the effect of viscosity on the flow in a bl.ide-to- 
blade surface of constant radius, which may be 
assumed to occupy a miu-span location since the 
effects of the hub and tip regions are neglected. 

The investigation has been confined to subsonic 
flows, but this limitation is not inherent in the 


method developed here. For the present discussion 
the solution surface can be considered to be the 
flat, two-dimensional region of a rectilinear cas- 
cade, with cambered blades of zero thickness. How- 
ever it will eventually be necessary to imagine this 
flat solution region as being wrapped onto the sur- 
face of a rotating right circular cylinder. Die 
introduction of viscosity into the calculation is 
accomplished by means of a vi: cous- inviscid inter- 
active calculation procedure. 

The inviscid calculation consists of a time- 
marching solution of the Fuler equations by MacCor- 
mack's Method. The viscous calculation proceeds in 
boundary layer and wake regions, and it solves a 
system of equations analogous to the set obtained 
by Horlock and Wordswortlv for the incompressible 
boundary layer on a helical blade. Although the 
viscous calculation is carried out on a cylindrical 
surface, the governing system of equations allows 
for the presence of a radial velocity component 
normal to that surface. The interaction between thi' 
viscous and inviscid calculations is accomplished 
by means of an iterative process. An iterative 
approach to the subsonic interaction problem is not 
unconnon; several researchers (c.g.. Ref. 17-18) 
have used this approach, in conjunction with the 
displacement thickness concept^" to obtain higlier 
approximations to flows. However, the present in- 
teractive method differs from these procedutes in 
two ways. First, the present method does not rely 
solely on the mechanism of a physical displacement 
of the outer flow streamlines by the viscous layer, 
to achieve coupling of the viscous and inviscid 
calculations. The interaction takes the forn. of an 
injection at blade surfaces (suction in the wake), 
but it is different from the usual source-sink dis- 
tribution technique in that this injection h..s a 
momentum and enthalpy character. St''ond, the appli- 
cation of boundary conditions to the viscous calcu- 
laion, and the viscous calculation itself are 
carried out in a manner suggested by the theory of 
matched asymptotic expansions. The details of the 
interactive procedure are discussed in Section IV. 

Tho viscous- invlsc id interactive calculation 
procedure which is described in this paper was used 
to calcu1.ite compressor flows for both rotor and 
stator passages. Some results of these niunerical 
calculations are presented. 

II . The Inviscid Solution 

For the inviscid solution we consider the invis- 
cid, rotational flow in a rectilinear cascade of 
zero thickness airfoils. The (a,0) coord iiiatf sv- 
stem used for this calculation (Fig. 1) is related 
to Cartesian coordinates by the relations, 

a “ X , P « y - B(x) (1) 

Also seen in Fig. 1 are the velocity components u 
and u , which are related to the Cartesian components 
of velocity in the following manner; 

u “ u cos O-uslnO (2) 

X s n 


u • u sin 6 + u cos 9 (3) 

y s n 

For this coordinate system the lime dependent 
Euler equations may be written in vector form as. 
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2»U . <iF . iC . „ . _ 


where , 


„ CU U + p CCS 0 
F » X 8 

pu u - p sin 3 
X n 

(e + p) u_ 


Ipu sec 0 
" 

Ipu u sec 0 

i V 

(,pu^ + p) sec 0 

> + p) u sec 6 
J n J 


0 ' 
-(pu^Up - p slnO]^ 

(PVs+Pcose)^ 


with e 


pLCyT + j (u^ + Uy)] and p “ pRT. 


A steady state solution of these equations In the 
cascade is obtained by a numer leal line inarching 
solution using MacCoriruick's Method . 

boundary Conditions 

As is often the case with time marching solu- 
tions of the Euler equations, the treatment of 
boundary conditions here consumes a disproportion- 
ately large part of the effort expended in the 
numerical solution. & careful treatment of certain 
boundary conditions in such problems seems to demand 
approaches which are somewhat involved. The dis- 
cussion of boundary conditions which appears in 
this section deals entirely with conditions as they 
exist for the first inviscid solution. Modifications 
to these boundary conditions are required for sub- 
sequent inviscid solutions within the viscous- 
inviscid iteration scheme, and a discussion of these 
modifications is deferred until Section IV. 

In a time marching solution of the Euler Equa- 
tions for the flew through a cascade (Fig. 1), 
essentially three different types of boundary regions 
are encountered. First, boundaries at which period- 
city conditions are the proper specification. 

This is the case at those portions of the boundary 
which connect the leading and trailing edges of the 
blades to the upstream and downstream boundaries. 
Second, the upstream and downstream boundaries which 
In this investigation are subsonic and permeable. 

It appears necessary that the specification of 
boundary conditions at these locations be compati- 
ble with the passage of wavelike disturbances 
through the boundary^O. Consequently, we treat the 
boundary conditions at these locations using the 
method of characteristics as suggested by Morettl^*. 
The details of this treatment are found inKeference 
22. Finally, solid wall boundaries represent a 
third type of situation. As the treatment of solid 
wall bouni’ary conditions is altered in subsequent 
inviscid solutions within the interactive scheme, 
we describe the situation during the first inviscid 
soiu. ion in scwie detail, so that the changes made 
for later solutions will be nxiro app.irent. 

Consider the numerical grid network iiear a 
blade surface which is depicted in Fig. 2. The 
grid lines J“2 ard J“3 He in the interior of the 
solution region, and solution values at these lo- 
cations are obtained fnxn the MacCormack algorithm. 
The J*1 grid line is a dummy point location, and it 
is at this location that the boundary ro.tuitlons 

ORK'JIXAI, P.\(, , i 
UP ()i; VM ! ' ’ 


are applied. The impermeability of the wall gives 
three relations, since three components of the 
vector G fsee Eq. (4)) are identically zero. Ihe 
effected components are Gi, G2 , and g 4 , where the 
numbers correspond to the position of the component 
within the vector. These relations are, 

“‘i.i • 

“l.t • -“l,! <‘> 

“l.l -“l,2 

The remaining component, G3 , reduces to (p sec 3). 

A fourth relation is obtained by again appealing 
to the method of characteristics. Fol lowing Moret 1 1 
, we seek to resolve those waves which propagate 
in a direction normal to tho boundary, apart from 
a translation tangent to the onundary due to the 
gross motion of the fluia. Using this approach we 
obtain (see Ref. 22) the compatibility relation, 

p - peUj, “-pc^c (u^ cos 0) - cu^ cos 0 ^ 
d0 

- u u 5-, (8) 

8 X d“' ' 

where c is the speed of sound, and the primes Je.notu 
differentiation in a direction defined by, 

^ ■ u + c sin 9, ^ “ sec 9(u - c) (9) 

dt X dt ' n 

Equation fS) may be integrated along a bicharactc- 
ristic line defined by Eq. (9) to obtain the wall 
pressure. Having obtained the wall pressure in 
this manner, a fourth relation is then available of 
the form, 

"l.l ■ • “l,2 

As a final note in this section, we mention 
that the Kutta condition is applied at the trailing 
edge of the blades by enforcing flow tangency. 


Extension to the Annular Cascade 

As a preliminary to the discussion of the vis- 
cous solution, it is worth noting that the numeri- 
cal solution for the inviscid flow in a rectiliunr 
cascade can be related to the flow in an annular 
cascade in a fairly simple way. To extend the 
results of the previous solution to the flow on a 
surface of constant radius which is in a state of 
radial equilibrium (l.e., zero radial velocity) and 
rotates about its axis, it is merely necessary to 
Imagine that the flat solution field is wrapped onto 
the surface of a rotating right circular cylinder. 


Ill . The Viscous Solution 

In this section, we develop the viscous equa- 
tions appropriate to the flow past cambered, yet 
strictly radial, blades. These blades and the co- 
ordinate system used in this development are seen 
in Fig. 3. The ((,C>^) coordinate system used in 
this section is shown in relation to a Cartesian 
coordinate system (X,Y, Z), and a cylindrical co- 
ordinate system (r,9,Z). The ^ and C coordinate 
lines are shown on a cylindrical surface fll“r" 
const.). The angle t((), which is measured on that 
surface, is the angle between the { coordinate line 


and a noncrator of tlu- aurface. If wo doacrlbo tho 
; coordinato lino aa a hollx-llko curve of anglo 
thon tho C cix)rdinato lino will bo a holix-liko 
curve of an«lo (270«’ + ,). rho * coordinate linos 
are atralKht lines normal to tho surface of tho 
cylinder. Finally, the velocity components In the 
and *1 directions are denoted by u,v and w, 
respectively. 

The curvilinear coordinates are related to Car* 
teslan coordinates In tho fol lowing way ; 

(«* 

Z ■ „ ^ cos > d5 + C sin V (ID 

X “ T) cos Y«7)sln(F (12,13 

where, , 

p sin i)d^ - C cos ji 
J o 

cp . „ 

Following Horlock and Wordsworth^, we confine 
our attention to the blade boundary layers which 
develop In a system that rotates about the Z axis 
with an angular velocity!;, and make some specifi- 
cation and assumptions. 

(1) Radial equilibrium Is specified for the 
extemal flow (w^ • 0). The ® subscript Indicates 
a location where C ’arge. 


+ - P„,) + sin t (po - p^u^) + y(r- 


oh ‘’h _ , . /Ou,^ . /i, 3T, 

pu ™ + pv — s “ u(p -^) +„(— ) + j., (k n) 


(17) 


The radial momentim equation (lb) Is uncoupled from 
the other three equations In the sense that the 
radial velocity (w) and the radial coordinate C,) 
do nut appear In Eqs . (14, lb, or 17). Consequently, 
if we confl.ie our attention to a surface of constant 
radius (Tj^cunst.), It Is possible to solve this 
system of equations with the appropriate boundary 
conditions by nuoierlcal marching techniques'^. To 
complete this system of equations, we take as the 
equation of state, 

p - pRT (18) 

and take the enthalpy, viscosity, and thermal con- 
ductivity r j be governed by the relations; 


T 

(19) 

P 

(T/T )’^^ 
ref' ref 

(20) 

(T/T )'®^ 
ref' ref' 

(21) 


Eqs. (14-21) are also the governing equations for 
the viscous wake. 


(11) The boundary layer thickness Is small 
compared to the blade chord; 

^ ^'(6., ) where 6 « 1. 

C N N 

(111) The chordwlse curvature of the blade Is 
of order (^). This Implies that 

C 5^ - (?(1), 
and that C ||r 

tlv) The chord Is small compared to the radius. 


It should be noted that, while the Invlscld and 
viscous calculations are both carried out on the 
same cylindrical surface, the coordinate systems 
used In these two calculations are different and the 
numerical grid systems would not In general coincide 
nor have the same orientation. 


IV. The Interactive Procedure 

The Interactive procedure takes the form of an 
Iteration between viscous and Invlscld solutions. 

In general teriw, this Iterative procedure Is as 
fol lows ; 


^ - (j) where e I . 

(v) The blade speed and u are of like order; 


~ AD • 


For a turbomachine It Is expected that <*■ * . 

With the ordering procedure established here. 

It la possible to reduce the equations of motion 
for the helical coordinate system to the appropriate 
boundary layer equations. The details of this re- 
duction may be found In Ref. 22. The boundary lay- 
er equations, correct to <f(t), which result are as 
follows : 

^ (pu) + ^ (pv) ■ 0 (14) 


(1) An Invlscld solution for the entire flow- 
field Is performed, with the appropriate boundary 
conditions . 

(11) Using boundary conditions, obtained from 
the Invlscld solution along blade surfaces and the 
wake centerline, the viscous calculation Is carried 
out. With the viscous ca Iculat Ion completed , certain 
adjustments are made In the Invlscld solution's 
boundary conditions, to reflect the presence of 
viscous layers. 

(ill) Steps (1) and (11) are repeated until an 
acceptable degree of convergence Is obtained. 

We now attend to the actual form of this Interaction. 


The Euler equations for steady flow nuiy be 
written. In Cartesian coordinates. In the vector 


form 


dif 0 ^ 


(22) 


du du 

5? P'' 5! 

dw dw 

pu + pv 5 ^ 


d du 

nr + 


/ 2 2 . , 2 ^ 
(pu - p„u ) sin 


(ID 


where 
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■pu^. 

2 ^ 

pu^ + P 

- 


PV'r 

2 ^ 

pu^ + P 

u (e + p) 


u (e + p) 

L 


1- / J 


( 16 ) 



Also, the ateady Navler-Stokea equatlona may be 
written In the vector form 


»r»? : ' 

- , du Oil 

p ^ T?' • ^ tJ: 


where ; 


f - 


,du ^ 


(23) 


dq_ du du 

I u,;e + p - A 

• Ou OU 

■ Y''^* tt' si' 




bu 


pvi ■ ‘‘(tJ ^ 

du 

X 


g - J pu,^ + p 


Lsv* •* 


du 

1 
7 

/du du 


du<^ 


-du du,.^ du^ 

•*■ X ' ■ ^ - 


/ du 


du 


- U . 




+ k 


U' 


J 


Now, we consider the flew In the linnedlate 
vicinity of the wall, where the viscosity and ther- 
mal conductivity are Important. We suppose that 
this viscous strip Is sufficiently thin (compared 
to the radius of curvature of the wall), and that 
the chordwlse extent of the region under considera- 
tion Is for the present sufficiently small, so that 
we are Justified In covering this region with a 
Cartesian coordinate system (see Fig. 4). An 
exact representation of the flow In this region Is 
given by a solution of the Navler-S tokes equations. 
Let f and g be the vectors constructed from this 
solution. Also, we suppose that some solution of 
the Euler equatlona will provide a close approxlaia- 
tlon of the exact solution when y > i, and let 3 
and jf be the vectors constructed from this Invlscld 
solution. Having Identified the vectors f, g, S’, 
and ^ with these two solutions, Eqs. (22) and (23) 
may be Integrated from y“0 to y“6, to give; 


^6 


J 


d • 

r» j 


3 dy 


(24) 


«6 ■ 


i 

5v 


f d4. 


(25) 


Since the two solutions are taken to coincide for 
<u> t, wa may specify ■ gj. Eqs. (24) and (25) 
■ay be combined then to give. 






(3 - f) dy 


(26) 


ORIGIN' /\L 1-. 

OF POOR UUALiil 


Eq. (26), which relates the two hypothetical solu- 
tions, will serve as a starting point for our dis- 
cussion of the solution technique In the vlsci>us 
layer. 

It Is not our Intention to solve the Navler- 

Stokes equations, therefore, we seek a suitable 

approximation of f and g on the Interval 0 ■ y t. 

We represent the exact solution by a composite 

function, f , where 
c 

f - 3 + f. (0) - 3 (27) 

C DO 

These functions are shosm In Fig. 4; 1^(0) corres- 
ponds to a boundary layer solution carried out using 
Invlscld values at <0* 0 as boundary conditions. 

The composite function f Is constructed In the 
spirit of a matched asymptotic expansion. The 
function f was chosen as an approximation of the 
exact solution for two reasons. First, we expect 
that this approach will have greater accuracy than 
the usual boundary layer solution. Second, we 
employ f^ because It has distinct computat lonn 1 
advantages within the context of the Iterative pro- 
cedure . 

Applying Eq. (27) to Eq. (26) gives, 

A ‘ 

/^o “ <*b(°»o s; j’ (=^o - “r 

o ' 


Eq. (28) can be used as the basis for an Iterative 
solution technique In the following way: 

(1) (/(D^. (/.'2)j, and (/f4) are Initially 

set equal to zero, and an Invlscld solution Is 
carried out. 

(11) Using Invlscld values at ^0, a boundary 
layer solution Is performed. 

(Ill) Using values obtained from the Invlscld 
and boundary layer solutions. Eq. (26) Is solved 
for new values of (^1) , and (/^) . The 

vector component 3) which contains the surface 
pressure Is evaluated using the method of characte- 
ristics. Since Eq. (28) requires only surface 
values from the Invlscld solutl, .i, a minimum of 
Interpolation Is required between the viscous and 
Invlscld grid systems. 

The Interaction model which has been described 
here can be conveniently used with those Invlscld 
solution procedures, which are currently employed 
to solve the Euler equations In primitive variable 
form. There Is an alternative method for dealing 
with the nusterlcal viscous- Invlsc Id Interaction 
when the Invlscld flow Is rotational, that being 
the displacement thickness approach, but It Is not 
conveniently used In a problem Involving complica- 
ted geometries. In such an approach, bodies are 
physically thickened, and It would be necessary to 
recompute the geometry of the problem at each step 
In the Iteration. In the present method, the 
geoewtry of the solid surface must be dealt with 
only once, and remains unchanged throughout the 
Iterative process. 

The form which the Interaction takes In the 
blade wake Is similar to the case of a wall boun- 
dary layer, which has been described. The details 
of the wake calculation are not reproduced here; 
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I ihoHc thf rcailfi' l» rcforrod to Ri-ffrrncc 22. 

V • NtmuMlcrtl RfsultH 

In thlM Mfi'ilon wi- prcHcnl himiu- lonult' ohtatn- 
i-J by applying thin Intorac 1 1 vf '’aKulat ton procf- 
tlurt' to thi' flow in a caHiado of zoro* th IcknotiH 
airfoilM, wlii<ii(' Hliapo 1 h that of a NACA • a ■ ,4 moan 
lino, with Cj j ■ .1 i»oo Rof. 24). TIu' MtaRRcr 
ant;l(' of tho caacado is I5‘’, Iho chordloiiRth of th»’ 
blades is i) ft., and the blade spacinR is .2 ft. 

We consider a subsonic fUmr of air throiiRli this 
cascade. The calculation lakes place on a cylindri- 
cal surface of radius 2.5 ft., and we consider both 
a rotor passaRe taiiRular velocity, 200 rad/sec), 
and a stator passage (zero angular velocity). Dim- 
to the nature of the equations which we are solving 
these two cases will differ >mly in the radial 
veliH itles in the boundary layers and wake witich 
result . 

In Figs. 5-8, values of the stre.,nMi8e velocity 
(it ) and density (p ) are plotted along several 8 ■ 
const, grid lines, where the location of the suction 
surface corrc'sponds to i * 0 ft., and the location 
of the pressure surface corresponds to 8 * .2 ft. 

The viscous- Invlsc Id Iterative scheme was run for 
four global iterations, and values from the first 
and last inviscid solutions are seen in these figu- 
res, The leading edge is denoted as I. .F., and the 
trailing edge by T.E. A noteworthy feature of these 
plots may be seen by comparing the first and last 
Inviscid solution values. Once the injection- 
suction boundary conditions are applied, an effec- 
tive bluntness is introduced at the leading edge. 
Also, an effective dlsplacemc-nt binlv surrounds the 
trailing edge. The new situation is nunu'rlcally 
less severe, and it may be' seen that a snu. 11 wavl- 
ness in the solution, which is apparent upstream 
and downstream of the blades, now disappears. Ihe 
final solution values at the leading edge behave 
as though a stagnation |>olnt had develo|H*d in the 
vicinity . 

An important advantage of an Interactive calcu- 
lation over a single inviscid calculation with a 
boundary layer added can he seen bv comparing the 
flow in the li.medlate vicinity of the trallingedge 
in Figs. 5 and b. At this location the interaction 
between the viscous flow and the inviscid flcW is 
strong, and the shape of the velocity profile 
changes slgnif ic,.nt ly between the first and last 
inviscid solutions. In Fig. b, a rapid dc>ce lerat ion 
of the fluid is indicated slightly dminstream of 
the trailing edge. A single Inviscid solutionwith 
a boundary layer added would not resolve this 
behavior. 

In Figs. 9 and 10, we plot rotor and stator 
radial velocity profiles for pressure surface 
boundary layer and wake locations. The plots in 
Figs. and 10 are taken from the final (fourth) 
viscous solution at locations about one third of a 
chord length behind the leading and trailing edge's 
respc'ct Ive ly . It may he seen that there is a large 
difference between the prtiflles obtained for a 
rotor passage and a stator passage, at both wake 
and boimdary layer locations. For a rotor passage 
the velocities are radially outward, and for a 
stator passage the velocities are radially inward. 
Also, it may he noted that generally larger values 
of the radial velocity are obtained in the wake 
than in the blade boundary layers. 


As an illustration of the cisnputer program's 
successful oj^ierat liin , values of the displacement 
thickness (6 ) are plotted over a portliin of the 
suction surface, for each of the four vlsciius solu- 
tions (Fig. II). We have limited thi chordwise 
extent of th)' ri'glon undej^ consideration in older 
to expand the vertical (6 ) scale, so that the 
convergence characteristics of the global Iterative 
scheme would he clearly visible. The abscissa in 
Fig. II corresponds to distance along the blad,' 
surface, nu'asurc'd from the leading edge. Ilie be- 
havior of successive solutions in Fig. II indicates 
convergence. Also, it appears that this conver- 
gence takes place quite rapidly, since the third 
and fourth solutions are virtually indistinguish- 
able even at this expanded scale. 

VI . Discuss Ion 

Hie numerical results of the preceding section 
Were taken from two solutions (rotor and stator) 
which Were carried ixit on an inviscid grid with 
90 X 20 disiensious. 'I1ie two ca Iculat ii>ns , which 
each required about 22 minutes (C.F.U. time) on a 
UNIVAC IIIO, were run for four global Iterations. 

The inviscid calculation procedure accounted for 
most of the run t It.ie . 

Hie computer program which has been developed 
in the court, of this study is currently limited 
in its ability to simulate real compressor fliws 
by the Idealizations which have been made. Ideali- 
zations such as blades of zero-thickness and 
strictly laminar fUiw have been intr>Hlu.ed to 
simplify the CkWiputat Iona 1 problem, but it is 
Important to note that these idealizations are not 
Inherent in our general approach to the viscous- 
Invisc id Interaction. The interactive calculation 
procedure which is presented here does not rely 
for its successful operation on the geonu'trical 
simplifications which have been made, and even 
depends very little on the precise form of the 
viscous and inviscid solutions. For example, an 
integral boundary layer calculation ctu Id be sub- 
stituted for the present viscixis marching procedure, 
or an implicit t lim> marching algorithm used to 
solve t h«' inviscid equations, and the overa 1 1 nature 
of the interactive calculation would not be much 
affected. This interactive scheme is novel in that 
it does not rely solely on the boundary layer dis- 
placement thickness, but incorporates Information 
related to momentum and enthalpy thicknesses as 
Well. The form of the interactive calculation 
procedure conveniently accomodates inviscid solu- 
tion procedures which are currently used to solve 
the F.uler equations in primitive variable form, 
and appears to have certain ctxnpu ta t iona 1 advantages 
for dealing with the viscous- invlsc id interaction 
when the inviscid fliw is rotational. 
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Figure 2. - Grid lines in the neighborhood 
of a blade surface. 
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Figure 3. ~ The coordinate system for the viscous solution. 
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Figure 4. - The inviscid, boundary layer, and 
composite solutions. 



Figure 5. - Chordwise velocity ifirst inviscid solution). 
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Figure 6. - Chordwise velocity (fourth inviscid solution). 
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Figure 8. - Density (^rth inviscid solution). 
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